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Abstract. When Newton's method, or Hahey's method is used to approximate the 
pth root of 1 — z, a sequence of rational functions is obtained. In this paper, a beautiful 
formula for these rational functions is proved in the square root case, using an interesting 
link to Chebyshev's polynomials. It allows the determination of the sign of the coefficients 
of the power series expansion of these rational functions. This answers positively the 
square root case of a proposed conjecture by Guo(2010). 



1. Introduction 

There are several articles and research papers that deal with the determination of the 
sign pattern of the coefficients of a power series expansions (see |5],[8] and the bibliography 
therein). In this work we consider this problem in a particular case. 

Let p be an integer greater than 1, and z any complex number. If we apply Newton's 
method to solve the equation = 1 — z starting from the initial value 1, we obtain the 
sequence of rational functions {Fk)k>o in the variable z defined by the following iteration 



Similarly, if we apply Halley's method to solve the equation = 1 — z starting from 
the same initial value 1, we get the sequence of rational functions {Gk)k>o in the variable 
z defined by the iteration 

It was shown in |3j and more explicitly stated in [4j that both Fk and Gk have power 
series expansions that are convergent in the neighbourhood of 2; = 0, and that these 
expansions start similar to the power series expansion of 2; 1— )■ {/I — z. More precisely, the 
first 2^ coefficients of the power series expansion of are identical to the corresponding 
coefficients in the power series expansion of z t— t- ^\ — 2;, and the same holds for the first 
3^ coefficients of the power series expansion of G^. It was conjectured [21 Conjecture 12], 
that the coefficients of the power series expansions of {Fk)k>2, {Gk)k>i and z 1— )■ ^1 — z 
have the same sign pattern. 
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In this article, we will consider only the case p = 2. In this case an unsuspected 
link to Chebyshev polynomials of the first and the second kind is discovered, it will allow 
us to find general formulae for and Gk as sums of partial fractions. This will allow 
us to prove Guo's conjecture in this particular case. Finally, we note that for p > 3 the 
conjecture remains open. 

Before proceeding to our results, let us fix some notation and definitions. We will 
use freely the properties of Chebyshev polynomials (T„)„>o and (f/n)„>o of the first and 
the second kind. In particular, they can be defined for a; > 1 by the folmute: 

\/ (f G M, T„(cosv9) = cos(n(^), and f/„ = — 'jp^n- (1) 

For these definitions and more on the properties of these polynomials we invite the 
reader to consult [B] or any general treatise on special functions for example Chapter 
22], $13.3-4] or [71 Part six], and the references therein. 

Let i7 = C \ [1, +oo), (this is the complex plane cut along the real numbers greater 
or equal to 1.) For z G ^2, we denote by v^l — z the square root of 1 — 2; with positive real 
part. We know that z 1— ?■ y/l — z is holomorphic in VL. Moreover, 

00 

\/zeWA), yr^=5^A^z"^, (2) 

m=0 



V^>0, A^^ = — ^1-— P^V (3) 



where 

v-m > n X = 

(1 - 2m)22™ \ m 

The standard result, is that ([2]) is true for z in the open unit disk, but the fact that 
for n > 1 we have A„ = /i„_i — /i„ < 0, where /i„ = 2~^'^(^^) ~ l/y^vm, proves the 
uniform convergence of the series ^ Am-^"* in the closed unit disk, and ([2]) follows by 
Abel's Theorem since z 1— t- \/1 — z can be continuously extended to -D(0, 1). 

Finally, we consider the sequences of rational functions (Ki)n>o, (-^n)n>o and {Gn)n>o 
defined by 

Vn^Az) = ^~'^Jr['\ V,{z) ^ 1, (4) 



^n+l(^) = ^FnW + ;^), Fo(^) = l, (5) 

Gn+i{z> = . , , , Go{z) = 1, (6) 



where (F„)n>o and (Gn)n>o are Newton's and Halley's iterations mentioned before, (in 
the case p = 2.) Since the sequence {Vn)n>o is simpler than the other two sequences, we 
will prove our main result for the V^'s, then we will deduce the corresponding properties 
for Fn and G„. 
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2. The Main Results 

We start this section by proving a simple property of that shows why it is sufficient 
to study the sequence of V^s to deduce the properties of Newton's and Halley's iterations 
the F„'s and G„'s : 

Lemma 2.1. The sequences (Ki)„, (-Fn)n {Gn)n of the rational functions defined 
inductively by and satisfy the following properties : 

(a) For z & Q and n > 0, we have : 



(b) For n > we have Fn = V21-1, and Gn = Vsn_i. 

Proof. First, let us suppose that 2; = x G (0, 1). In this case, we see by induction that 
all the terms of the sequence (Ki(a;))„>o are well defined and positive, and we have the 
following recurrence relation : 



Vn+l{x) - yi -X ^ 1-X + Vnjx) - y/l - x{l + Vnjx)) 
Vn+l{x) + y/1 - X 1-X + Vn{x) + y/1 - x{l + Vn{x)) 

1 — \/l — X Vn{x) — ^/l — X 
l + ^/l-X Vn{x) + ^/l-x' 

Therefore, using simple induction, we obtain 

V.„>0. ^"W-^=(^-^V". (7) 

Now, for a given n > 0, let us define R{z) and L{z) by the formulae : 



Clearly, R is analytic in f2 since — z) > for z G ^l. On the other hand, if 

Vn = An/ Bn whcrc An and Bn are two co-prime polynomials then 



{An{z) - VT^Bn{z)f 

AUz)-il-z)Bl{z) ■ 

So L is meromorphic with, (at the most,) a finite number of poles in Q. Using analyticity, 
we conclude from ([7]) that L{z) = R{z) for z E fl\V for some finite set V C fl, but this 
implies that the points in V are removable singularities, and that L is holomorphic and 
identical to R in Q. This concludes the proof of ([a]). 
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To prove (jb]), consider 2; = a; G (0, 1), as before, we see by induction that all the 
terms of the sequences {Fn{x))n>o and (Gn(a^))n>o are well-defined and positive. Also, we 
check easily, using and (E]), that the following recurrence relations hold : 

Fn+i{x) - y/1 - X _ / l-Vl-a: \^ Fnjx) - yi - X 
Fn+i{x) + ~ \l + ) ' F^{x) + VT^' 

and 



Gn+l{x) + + Gnix) + ^1 - X ' 

It follows that 

on 

Fnjx) - _ / I - VT^ \ _ V2«-l{x) - y/T^ 



Gn{x) — \/l — X _ / 1 — Vl — x\ _ Vz^-i{x) — ^Jl — X 
Gn{x) + ^l-X~ \l + ^l-x) ~ Vy.^i{x) + ^l-X 

Hence = V2"-i(x) and Gn{x) = V3n_i(x) for every x G (0,1), and (jb]) follows by 

analyticity. □ 



Lemma [2.11 provides a simple proof of the following known result. 

Corollary 2.2. The sequences {Vn)n, {Fn)n o^nd {G„)n of the rational functions defined 
inductively by (jl]), (jS]) and converge uniformly on compact subsets ofQ to the function 

Z I—)- a/1 — z. 



Proof. Indeed, this follows from Lemma \2A] and the fact that for every non-empty compact 
set C we have : 



sup 



1 + 



< 1, 



which is easy to prove. 



□ 



One can also use Lemma 12.11 to study Vn in the neighbourhood of 2; = : 
Corollary 2.3. For every n > 0, and for z in the neighbourhood of we have 



and 



where Am is defined by 



v;(^) = v/r^ + o(^"+i), 

n 

K(^) = 5ZAmz'" + 0(z"+i). 



(9) 



m=0 
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Proof. Indeed, using Lemma 12.11 we have : 



Vniz) = ■ 



[1 + v/r^)"+l + (1 - 



z 



in+1 



;i + vr^)"+i - (1 - vr^)"+^ 



(1 + v^r^)"+i - (1 - v/r^)"+i (1 + v/r^)"+ 



(1 + v/r^)2{«+i) - 



In particular, for z in the neighbourhood of 0, we have Vn{z) = a/1 — z + 0(z"^^), which 

is dHD. 

On the other hand, using ([2]), we obtain 



Vn{z) = J2^n.z'^ + 0{z''' 



m=0 

which is (El). □ 



Recall that we are interested in the sign pattern of the coefficients of the power series 
expansion of F„ and Gn, in the neighbourhood of 2; = 0. Lemma [2.11 reduces the problem 
to finding sign pattern of the coefficients of the power series expansion of Vn- But, Vn is 
rational function and a partial fraction decomposition would be helpful. The next theorem 
is our main result : 



Theorem 2.4. Let n be a positive integer, and let Vn be the rational function defined by 
the recursion (j4]) . Then the partial fraction decomposition of Vn is as follows : 

V{z) = l-'-, 

Vn{z) = 1 - - - — — — - V " Y ^fc N ' ^^"^ n>2. 

2 2{n + l) ^ l_2;cos2(^) 



Proof. Let us recall the fact that for x > 1 Chebyshev polynomials of the first and the 
second kind satisfy the following identities : 
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Thus, for < a; < 1 we have 



X 



l\ (1 + VI - + (1 - 



X 



2\n 



l\ (1 + Vl - - (1 + VT^^Y+^ 



xj 2^/1 -X' 

So, using by Lemma [2.1( a). we find that 

T/ (rr'i\ — x''~^^Tn+l{l/x) _ Pnjx) , , 

where P„(x) = x"+iT„+i(l/x) and g„(x) = x"f/„(l/x). 

It is known that f/„ has n simple zeros, namely {cos(fc^„) : 1 < k < n} where On = 
n/{n + 1). So, if we define 

A2™ = {l,2,...,2m} 

and 

A2„,-i = {l,2,...,2m-l}\{m}, 

then the singular points of the rational function Pn/Qn are {A^ : k G A^} with = 
sec{k6n). Moreover, since 

= A^+1 cos(7rA;) = (-1)'=A^+V 0, 
we conclude that these singular points are, in fact, simple poles with residues given by 

Rp« /^^" X ^ - x3 ^n+l(l/Afc) 

But, from the identity Unicosip) = sin((ri + simp we conclude that 

cos(^sin((n + cos{{n + l)<p) 
Unicos^) = — (n + 1) — 



hence. 



and finally. 



sm (f sm ip 



U:^{cos{ken)) = {n + l) 



sm^ikOnY 



Res ( ^,Xk] = ^^sm\k9n) = ^^tan^ikOn). (13) 



Qn J n + 1 n + 1 
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From this we conclude that the rational function i?„ defined by 



RJx) 



+ 1 ^ 



Qnix) n + 1 X- Xk 

Pn{x) _^ 1 tan2(A;6'„ 



+ 1 i-xcoM„y ^ ' 



Qn{x) n + 1 ^ 1 - xcos{k6n)' 

is, in fact, a polynomial, and degi?„ = degP„ — degQ„ < (n + 1) — (n — 1) = 2. 
Noting that k^-^n + l — kisa permutation of A„ we conclude that 

Etan^(A;^^„) tan^(A;^„) 
l-xcos(k9n) ~ ^ l + xcos(k9n) 

tan^(A;6'„) tan^(A;6'r, 



2 \ 1 + X cos(A;^^„) l — xcos(k6n) 

keAn 



and using the fact that 



- — x'^ cos^ (kOn) 



2 4 

2 2; COS 99 

1 + Z COS V9 + 



we find 



1 — Z COS^ 1 — 2; COS^ (p 



E 1 , = E tan2(fc^„) V sin2(A:^„] 

cos^(/c6'„) sin^(/c6'„) 



x'^ cos^ (kOn) 

k=l 

where we added a "zero" term to the last sum for odd n. 

Combining this conclusion with ( IT^ we conclude that there exists a polynomial Sn with 
deg 5"^ < 2 such that 

Pnjx) x^ -A sin^(2fc6'n) 

"^""^ " + 4(n + 1) ^ 1 - x2 cos2(A;^„) ' 

Moreover, S'„ is even, since both P„ and Q„ are even. Thus, going back to I IT2\i we conclude 
that there exists two constants a„ and /3n such that 

T./ N o sin^(2/c6'„) 

Kw=«„+A.+j(;^i:^— (15) 

But, from Corollary I2.3[ we also have Vn{z) = 1 — \z + Oi^z"^) for n > 1, thus a„ = 1 and 

= -1/2. 

Finally, noting that the terms corresponding to k and n + 1 — k in the sum f[T^ are 
identical, we arrive to the desired formula. This concludes the proof of Theorem 12.41 □ 
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Theorem 12.41 allows us to obtain a precise information about the power series expan- 
sion of Vn in the neighbourhood of z = : 

Corollary 2.5. Let n be a positive integer greater than 1, and let Vn be the ratio- 
nal function defined by (j4]). Then the radius of convergence of power series expansion 
Ylm=o^'^^^^'^ of in the neighbourhood of z = is sec^(^^) > 1, and the coefficients 
{Am^)m>o satisfy the following properties : 

(a) For < m < n we have Am^ = X^, where Xm is defined by 

(b) For m > n we have Am^ < and 

\ m ) V / n + l 

m=n+l ^ ' 



Moreover, for every n > and every z in the closed unit disk -D(0, 1) we have 

\Vniz) - VT^zl < ^ 





2 , 


< - 


\z\ 




\ Tin 



In particular, (Ki)n>o converges uniformly on D{0, 1) to the function z h- )■ \/l ~ z. 



Proof. Let us denote Tr/{n + 1) by By Theorem 12.41 the poles of Vn, for n > 1, are 
{sec^(fc6'„) : 1 < k < \n/2\} and the nearest one to is sec^(6'„). This proves the 
statement about the radius of convergence. 



Also, we have seen in Corollary 12. 3[ that in the neighbourhood of 2; = we have 

n 

This proves that for < m < we have Am^ = Xm which is 

To prove ()b]), we note that for 1 < /c < [n/2\, and z G -D(0, sec^(6'„), we have 



L^ = ^cos2'"(A:^„ 

1 — Z COS^ kUn 



m=0 

hence 



00 / [n/2i 

Vn{z) = 1 - I - E ^7;^ E cos^'-lfc^n) sin2(2A:0„) 

?Ti=0 \ ^==1 



This gives the following alternative formulae for Am when m > 2 : 

^ L"/2J 

AL") = cos2(™-i)(A;^„)sin2(A;^„), 



k=i 



^f^cos2(™-i)(A;^„)sin2(fc^„), (16) 



n + _ 

k=i 
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where ( IT6|) is also valid for m = 1. This proves in particular that Am'^ < for m > n. 

dm=0 



Since the radius of convergence of Y^^=o^rnz'^ is greater than 1 we conclude that 



(17) 



m=0 



but, we can prove by induction from (j4]) that = l/(n + 1). Therefore, (1T7|) implies 

that, for n > 1 we have 

oo ^ n n 



m=n+l 



m=0 



m=l 



1 " 

n + 1 



n + 1' 



where /x^ = 2-2"^ (2™), which is ©. 
For n > 1 and z G -D(0, 1) we have 



m=0 



< 



E 4 



n) m 
m 



m=n+l 



< k 



< 



E (-4 



n)\ 
m ) 



m=n+l 



1 f2n 



)2n 



n + 1 



in+l 



On the other hand, using ([2]) we see that for n > 1 and z G -D(0, 1) we have 



m=0 



< 



m=n+l 



< k 



n+l 



E ^/"'^ 



-1 ~ ^^r, 



m=n+l 



< 



1 /2n 



)2n 



in+l 



Combining ([TK|) and ffT^ . and noting that 2 ") < l/v^vm we obtain 
which is the desired conclusion. 



:i8) 



(19) 



□ 



The following corollary is an immediate consequence of Lemma 12.11 and Corollary 
12.51 It proves that Conjecture 12 in [3] is correct in the case of square roots. 
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Corollary 2.6. The following properties of {Fn{z))n>o and {Gn{z))n>o, the Newton's and 
Halley's approximants of a/1 — z, defined by the recurrences and (jH]) hold : 

(a) For n > 1, the rational function Fn{z) has a power series expansion 1 — Ylm=i Bm'^ z"^ 
with sec^(2~"7r) as radius of convergence, and > for every m > 1. Moreover, 

2 '-'^ 



VzG 1^(0,1), \Fn{z) - VT^\ < ^ 



V2F~^' 

(b) For n > 1, the rational function Gn{z) has a power series expansion 1 — X]m=i Grn^ z" 
with sec^(3~"'7r) as radius of convergence, and Ct^ > for every m > 1. Moreover, 



\fzeD{0,l), \Gn{z) - VT^l < ^ ■ 



\z\ 



TT - 1 

Remark. Guo's conjecture Conjecture 12] in the case of square roots is just the fact 
that Bm^ > for n > 1, m > 1, and that Cm^ > for n, m > 1. The case of pth roots 
for p > 3 remains open. 

Conclusion. In this paper, we presented a link between the rational functions approxi- 
mating z I—)- a/1 — z obtained from the application of Newton's and Halley's method, and 
Chebyshev Polynomials. This was used to find the partial fraction decomposition of this 
rational functions, and the sign pattern of the coefficients of their power series expansions 
was obtained. Finally, this was used to prove the square root case of Conjecture 12 in [3], 
and to give estimates of the approximation error for z in the closed unit disk. 
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